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ABSTRACT : The design of composite structures requires an evaluation of their safety 
and durability under service loads and possible overload conditions. This paper presents 
a computational tool that has been developed to examine the response of stiffened compos- 
ite panels via the simulation of damage initiation, growth, accumulation, progression, and 
propagation to structural fracture or collapse. The structural durability of a composite 
panel with a discontinuous stiffener is investigated under compressive loading induced by 
the gradual displacement of an end support. Results indicate damage initiation and pro- 
gression to have significant effects on structural behavior under loading. Utilization of an 
integrated computer code for structural durability assessment is demonstrated. 

KEY WORDS: composite materials, composites, computational simulation, damage, 
degradation, durability, fracture, laminates, simulation, stiffened panel, structural 
degradation. 


INTRODUCTION 

S tiffened COMPOSITE PANELS are used in many structural components to sat- 
isfy requirements of reduced weight, increased stiffness, and stability. Design 
considerations with regard to the durability of stiffened panels require an a priori 
evaluation of damage initiation and propagation mechanisms under expected 
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service loading. Concerns for safety and survivability of critical components re- 
quire a quantification of the composite structural damage tolerance during over- 
loads. 

The numerous options available in the tailoring of composite structures make 
composites more versatile for fulfilling structural design requirements. However, 
these same design flexibilities render the assessment of composite structural re- 
sponse and durability more complex, prolonging the design and certification pro- 
cess and adding to the cost of the final product. It is difficult to evaluate compos- 
ite structures because of the complexities in predicting their overall congruity and 
performance, especially when structural degradation and damage propagation 
take place. The predictions of damage initiation, damage growth, and propaga- 
tion to fracture are important in evaluating the load carrying capacity, damage 
tolerance, safety, and reliability of composite structures. Quantification of the 
structural fracture resistance is also fundamental for evaluating the durability/life 
of composite structures. The most effective way to obtain this quantification is 
through integrated computer codes that couple composite mechanics with struc- 
tural analysis and damage progression modelling. The Composite Durability 
STRuctural ANalysis (CODSTRAN) computer code [1J has been developed for 
this purpose. The simulation of progressive fracture by CODSTRAN has been 
validated to be in reasonable agreement with experimental data from tensile tests 
[2]. Recent additions to CODSTRAN have enabled investigation of the effects of 
composite degradation on structural response [3], composite damage induced by 
dynamic loading [4], composite structures global fracture toughness [5], effect of 
the hygrothermal environment on durability [6], structural damage and fracture 
in composite thin shells [7], structural durability of a composite pressure vessel 
[8], overall evaluation of damage progression in composites [9], damage and 
fracture of stiffened shell panels [10], and design implications of progressive frac- 
ture in composite shell structures [11]. The present paper demonstrates the capa- 
bility of CODSTRAN to evaluate the response and degradation of a discon- 
tinuous^ stiffened composite panel under a displacement controlled loading. 

THE CODSTRAN METHODOLOGY 

CODSTRAN is an open-ended computer code integrating selected modules on 
composite mechanics, damage progression modelling, and finite element analy- 
sis. The damage progression module [1] keeps a detailed account of composite 
degradation for the entire structure and also acts as the master executive module 
that directs the composite mechanics module [12] to perform micromechanics, 
macromechanics, laminate analysis and synthesis functions. It also calls the finite 
element analysis module [13] with anisotropic thick shell analysis capability to 
model laminated composites for global structural response. A convenient feature 
of the utilized finite element module is that structural properties are input and 
generalized stress resultants are output at the nodes rather than for the elements. 
The anisotropic generalized stress-strain relationships for each node are revised 
according to the composite damage evaluated after each finite element analysis. 
Subsequent to damage and composite degradation, the model is automatically 
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Figure 1 . CODSTRAN simulation cycle . 


updated with a new finite element mesh and properties before the structure is 
reanalyzed for further deformation and damage. 

Figure 1 shows a schematic of the computational simulation cycle in 
CODSTRAN. The ICAN composite mechanics module is called before and after 
each finite element analysis. Prior to each finite element analysis, the ICAN 
module computes the composite properties from the fiber and matrix constituent 
characteristics and the composite layup. The laminate properties may be different 
at each node. The finite element analysis module accepts the composite proper- 
ties that are computed by the ICAN module and performs the analysis at each 
load increment. After an incremental finite element analysis, the computed gen- 
eralized nodal force resultants and deformations are supplied to the ICAN 
module that evaluates the nature and amount of local damage, if any, in the plies 
of the composite laminate. Individual ply failure modes monitored by 
CODSTRAN include the failure criteria associated with the negative and positive 
limits of the six ply-stress components (a m , <* 02 , cr /23 , tf/is) and a 

combined stress or modified distortion energy (MDE) failure criterion, and in- 
terply delamination due to relative rotation (RR) of the plies [12]. 

For the purpose of the present discussion, the following terminology is used to 
describe the various stages of degradation in the composite structure: (1) damage 
initiation refers to the start of damage induced by loading that the composite 
structure is designed to carry; (2) damage growth is the progression of damage 
from the location of damage initiation to adjacent regions; (3) damage accumula- 
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tion is the increase in the amount of damage in the damaged region with addi- 
tional damage modes becoming active; (4) nodal fracture is a through-the- 
thickness fracture at a node because of major tensile or compressive failures in 
all plies of the laminate. Nodal fracture is usually a characteristic in precipitating 
rapid damage propagation toward structural fracture. 

STIFFENED COMPOSITE PANEL 

The demonstration example for this paper consists of a stiffened composite 
panel, depicted in Figures 2a and 2b, subjected to axial compression. The finite 
element model shown in Figure 2a uses thick shell elements with duplicate nodes 
where there are sudden changes in composite properties. The use of duplicate 
nodes to achieve accurate structural representation at locations of abrupt change 
in the laminate was discussed in a previous paper [10]. The composite system for 
this example is made of AS-4 graphite fibers in a high-modulus, high-strength 
epoxy matrix (AS-4/HMHS). The corresponding fiber and matrix properties are 
obtained from a databank of composite constituent material properties resident in 
CODSTRAN [12]. 

The properties used for the AS-4 graphite fibers and HMHS epoxy matrix are 
given as follows: 

AS-4 Fiber Properties: 

Number of fibers per end = 10000 
Fiber diameter = 0.00762 mm (0.300E-3 in) 

Fiber Density = 4.04E-7 Kg/m 3 )(0.063 lb/in 3 ) 



Figure 2a. Stiffened composite panel finite element model AS-4/HMHS [[0/± 45/90] ,]e. 
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Figure 2b. Stiffened composite panel AS-4/HMHS [[0/±45/90] s ]e cross section and plan 
(all dimensions are in mm). 


Longitudinal normal modulus = 227 GPa (32.90E + 6 psi) 

Transverse normal modulus = 13.7 GPa (1.99E + 6 psi) 

Poisson's ratio (y, 2 ) = 0.20 
Poisson's ratio (j> 23 ) = 0.25 
Shear modulus (G, 2 ) = 13.8 GPa (2.00E + 6 psi) 

Shear modulus (G 23 ) = 6.90 GPa (1.00E + 6 psi) 

Longitudinal thermal expansion coefficient = — 1.0E — 6/°C (— 0.55E-6/°F) 
Transverse thermal expansion coefficient = 1.0E — 5/°C ( + 0.55E — 5/°F) 
Longitudinal heat conductivity = 43.4 J— m/hr/m 2 /°C (580BTU — 
in/hr/in 2 /°F) 

Transverse heat conductivity = 4.34 J-m/hr/m 2 /°C (58 BTU — 
in/hr/in 2 /°F) 

Heat capacity = 712 J/Kg/°C (0.17 BTU/lb/°F) 

Tensile strength = 3,723 MPa (540 ksi) 

Compressive strength = 3,351 MPa (486 ksi) 

HMHS Matrix Properties: 

Matrix density = 3.40E — 7 Kg/m 3 (0.0457 lb/in 3 ) 

Normal modulus = 4.27 GPa (620 ksi) 
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Poissons ratio = 0.34 

Coefficient of thermal expansion = 0.72/°C (0.4E — 4/°F) 

Heat conductivity = 1.25 BTU -in/hr/in 2 /°F 
Heat capacity = 0.25 BTU/Ib/°F 
Tensile strength = 84.8 MPa (12.3 ksi) 

Compressive strength = 423 MPa (61.3 ksi) 

Shear strength = 148 MPa (21.4 ksi) 

Allowable tensile strain = 0.02 
Allowable compressive strain = 005 
Allowable shear strain = 0.04 
Allowable torsional strain = 0.04 

Void conductivity = 16.8 J-m/hr/m 2 /°C (0.225 BTU -in/hr/in 2 / °F) 

Glass transition temperature = 216°C (420°F) 

The HMHS matrix properties are representative of the 3501-6 resin. The skin 
laminate consists of forty-eight 0.132 mm (0.00521 in) plies resulting in a com- 
posite thickness of 6.35 mm (0.25 in). The fiber volume ratio is 0.60. The lami- 
nate configuration in [[0/ ±45/90], ] 6 . The 0° plies are in the axial direction of 
the stiffener, along the x axis indicated in Figure 2a. The width of the stiffened 
panel is 991 mm (39 in) and it has a length of 559 mm (22 in). The stiffener hat 
sections are made from the same AS-4/HMHS composite structure as the outer 
skin. The stiffeners are perfectly bonded to the skin at all surfaces of contact. 
This example is physically very similar to the stiffened composite panel analyzed 
in Reference [14]. Similar to Reference [14], large displacements are taken into 
account. The novelty in the present paper is that composite constituent level pro- 
gressive damage simulation is integrated into the structural analysis of the 
stiffened panel. 

Axially compressive loading is applied by imposing a gradually increasing uni- 
form axial displacement at the clamped edge of the panel. A displacement con- 
trolled CODSTRAN simulation is employed to monitor damage initiation and 
progression as the panel is loaded. Figure 3 shows the relationship between in- 
duced loading due to the imposed displacement of the clamped edge and the pro- 
duced damage during damage initiation and growth stages. First ply failure is in 
the surface ply of the skin near the end of the stiffener. 

The scalar damage variable, shown on the ordinate of the graph in Figure 3, is 
derived from the total volume of the composite material affected by the various 
damage mechanisms [5,7]. Computation of the shown scalar damage variable has 
no interactive feedback on the detailed simulation of composite degradation. 
CODSTRAN is able to simulate varied and complex composite damage 
mechanisms via evaluation of the individual ply failure modes and associated 
degradation of laminate properties. In general, the type of damage growth and the 
sequence of damage progression depend on the composite structure, loading, 
material properties, residual stresses induced by the curing process, and 
hygrothermal conditions. 

With reference to Figure 3, damage initiation is under a 635 kN (143 kip) load- 
ing by compressive failures in the 0° first ply at approximately 6.4 mm from the 
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Compressive Load (kN) 

578 623 667 712 756 



Compressive Load (kips) 

Figure 3. Damage initiation and growth under compressive load AS-4/HMHS 
[[0/± 4 5/90] s ] e- 


stiffener web ends toward the clamped edge of the panel. As loading increases, 
damage grows contiguously in the plane of the skin. Across the laminate thick- 
ness, damage growth is first into the zero degree plies. When loading reaches 737 
kN, the ±45° angle plies also begin to participate in damage progression. After 
the 737 kN loading level, damage growth is accompanied with negligible increase 
in loading. After a 737.4 kN ultimate load, corresponding to a cumulative damage 
level of approximately 0.14 percent, the induced compressive load decreases as a 
through-the-thickness damage propagation stage is entered with increasing global 
strain. The damage propagation mode is highly localized to the unstiffened com- 
posite skin areas immediately adjacent to the ends of the stiffener toe elements. 

Figures 4a and 4b show the Ply 1 longitudinal stress contours under a 689 kN 
load during damage growth. Composite properties are degraded according to the 
accumulated damage, resulting in the lowering of ply stresses in the damaged re- 
gions. Damage is localized to the skin, adjacent to the ends of the stringer toe ele- 
ments bonded to the skin. Due to the localized nature of damage, the overall pat- 
tern of stress contours shown in Figure 4a is not significantly affected during the 
initial stages of damage. The significant local effects at the damaged region are in- 
dicated in Figure 4b by the detailed stress relief contours. 

Figure 5 shows the overall relationship between the applied global strain and 
the corresponding total axial compressive force in the panel. The CODSTRAN 
simulation indicated in Figure 5, takes into account geometric nonlinearities 
caused by large displacements as well as the effects of structural degradation with 
damage. Linear analysis and geometrically nonlinear analysis results, not taking 




92 


Levon Minnetvan et al. 



Figure 4a. Ply 1 global longitudinal stress contours after initial damage growth AS-4/HMHS 
[10/ ± 45m], ]e 



Figure 4b. Ply 1 stress relief details in critical region after initial damage growth AS-4/HMHS 
[[0/±45m],] 6 . ' ' 
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Figure 5. Induced compressive load with end displacement AS-4/HMHS [[0/±45/90] s ]6. 


into account the effects of structural damage, are also shown in Figure 5 for com- 
parison with CODSTRAN simulation. Before the damage initiation stage, 
CODSTRAN simulation is identical to geometrically nonlinear structural analy- 
sis. However, the information provided by CODSTRAN with regard to damage 
progression and damage tolerance during the various degradation stages cannot 
be obtained from any traditional linear or nonlinear structural analysis result. 

Figure 6 shows the direct relationship between the applied global strain and the 
volume of damage produced in the composite structure. After the damage initia- 
tion stage, overall damage volume increases precipitously yet steadily with the 
applied global strain through the maximum loading stage and beyond. 

A measure of composite structural damage tolerance is obtained by defining the 
Strain Energy Release Rate (SERR) as the amount of incremental work done on 
the structure per unit volume of the created damage during degradation [5,7]. 
Figure 7 shows the SERR as a function of the applied global strain on the panel. 
The peak SERR level occurs at damage initiation, corresponding to a global 
deformation of approximately 4000 microstrains. SERR falls to a much lower 
level after damage initiation, indicating the lack of damage tolerance of this struc- 
ture after initial damage. Fluctuations in the SERR indicate changes in the struc- 
tural resistance against damage under progressive end displacement. The maxi- 
mum loading level of 737.4 kN (166 kips), corresponding to an applied global 
deformation of approximately 5\0Qpe, is marked by a small but sharp local peak 
in the SERR. 

Figure 8 shows the z components of displacements under a 4500pe global strain 
during damage growth. The overall displacement contours depicted in Figure 8 
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Figure 6. Damage propagation with end displacement AS-4/HMHS [[0/ ±45/90],],. 
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Figure 7. Strain energy release rate with end displacement AS-4/HMHS ([0/± 45/90] S J 6 . 
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DISPLACEMENTS AFTER DAMAGE INITIATfON 

Figure 8. Displacements (z components) after damage initiation AS-4/HMHS 
[[0/±45/90] a J e . 


are not affected by the localized damage. Figure 8 indicates that the apparent 
global response of the stiffened panel, as may be observed during a physical test, 
does not show any significant signs of the damage progression that is taking place 
in the plies of the laminate. 

Figure 9 shows the x components of displacements under the same 4500/*e 
global strain as in Figure 8. Axial displacement coutours shown in Figure 9 are 
only slightly affected near the local damage. Figure 9 indicates that in order to 
detect damage growth in a laboratory test it is necessary to measure accurately 
the local strains in the damaged zone. 

Figure 10 shows Ply 1 longitudinal stress contours immediately after through- 
the-thickness laminate fractures occur under a 737.4 kN (166 kip) ultimate load- 
ing. The stress relief contours near the stringer toe edges shown in Figure 10 are 
indicative of laminate failure locations. 

CONCLUDING REMARKS 

The significant results and overtones from this investigation in which 
CODSTRAN (Composite Durability STRuctural ANalysis) is used to evaluate 
structural response of a stiffened composite panel, considering damage initiation 
and progression effects, are as follows: 

1. Computational simulation, with the use of established composite mechanics 
and finite element modules, can be used to predict the damage tolerance. 
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Figure 9. Axial displacements after damage initiation AS-4/HMHS [[0/±45/90J,]e- 
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safety, and durability of built-up composite structures such as composite skin 
panels with integrated stiffeners. 

2. CODSTRAN adequately tracks the damage initiation, growth, and subse- 
quent propagation to fracture for any composite structure. 

3. Non-destructive evaluation of in-service structural integrity is facilitated by 
the prediction of damage initiation/progression locations and mechanisms. 

4. CODSTRAN simulations may be carried out under either load controlled or 
displacement controlled conditions. 

5. The demonstrated procedure is flexible and applicable to all types of constitu- 
ent materials, structural geometry, and loading. Hybrid composite structures, 
composite laminates containing homogeneous materials such as metallic foils, 
as well as binary composites can be simulated. The hygrothermal environ- 
ment, residual stresses induced by the curing process, local defects due to fab- 
rication error and/or accidental damage may be included in a CODSTRAN 
simulation. 

6. Fracture toughness parameters such as the structural fracture load and the 
ultimate load are identifiable for any structure by the demonstrated method. 
It is also useful to carry out CODSTRAN simulations prior to physical testing 
to guide the data acquisition strategy and to enable the detailed interpretation 
of experimental results with regard to damage initiation/progression 
mechanisms and damage tolerance. 

7. The availability of CODSTRAN facilitates composite structural design and 
certification by allowing the efficient and effective evaluation of design options 
and by providing early design loads. 

8. Computational simulation by CODSTRAN represents a new global approach 
to structural integrity assessment. 
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